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Abstract 



A version of the virial theorem is derived in a brane- world scenario in the framework of a warped 
DGP model where the action on the brane is an arbitrary function of the Ricci scalar, C{R). The 
' (— I '\ extra terms in the modified Einstein equations generate an equivalent mass term (geometrical 

^ j mass), which give an effective contribution to the gravitational energy and offer viable explanation 

to account for the virial mass discrepancy in clusters of galaxies. We also obtain the radial velocity 
dispersion of galaxy clusters and show that it is compatible with the radial velocity dispersion 



profile of such clusters. Finally, we compare the result of the model with C{R) gravity theories. 

PACS numbers: 04.50.-h, 04.20.Jb, 04.20.Cv, 95.35. +d 
^ \ Key words: Virial mass, DGP brane-world, Boltzmann equation, Radial velocity dispersion. 
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1 Introduction 

The intriguing possibility that our universe is only part of a higher dimensional space-time has raised 
a lot of interest in the physics community [1], motivated by the developments in superstring and 
M-theory. According to the brane-world scenario, all matter and gauge interactions reside on the 
brane, while gravity can propagate in the whole five-dimensional space time (figure 1) [2]. Higher 
dimensional models have a long history, but it revived by works of Randall and Sundrum (RS) in 
1999 [3]. They have introduced two models in order to solve the hierarchy problem in particle physics, 
however after a while these two models, because of their interesting properties, could attract a salient 
attention in cosmology. In their first model they suppose two brane which our brane has a negative 
tension (Shiromizu et al. [4] have shown that this model is unphysical). In their second model they 
suppose a brane with infinite extra dimension. In that model our universe has a positive tension. 
The effective four-dimensional gravity in the brane is modified by extra dimension. The cosmological 
evolution of such a brane universe has been extensively investigated and effects such as a quadratic 
density term in the Friedmann equations have been found [5, 6, 7]. An alternative scenario was 
subsequently proposed by Dvali, Gabadadze and Porrati (DGP) [8]. The DGP proposal rests on 
the key assumption of the presence of a four-dimensional Ricci scalar in the brane action. There 
are two main reasons that make this model phenomenologically appealing. First, it predicts that 
four-dimensional Newtonian gravity on a brane-world is regained at distances shorter than a given 
crossover scale (high energy limit), whereas five-dimensional effects become manifest above that 
scale (low energy limit) [9]. Second, the model can explain late-time acceleration without having to 
invoke a cosmological constant or quintessential matter [10]. An extension of the DGP brane-world 
scenario have been constructed by Maeda, Mizuno and Torii which is the combination of the RS II 
model and DGP model [11]. In this combination, an induced curvature term appears on the brane in 
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Figure 1: The confinement of matter to the brane, while gravity propagates in the bulk. 

the RS II model which has been cahed the warped DGP brane-world in hterature [12]. In this paper, 
we consider the effective gravitational field equations within the context of the warped DGP brane- 
world model where the action on the brane is an arbitrary function of the Ricci scalar, C{R), (which 
is called the modified DGP brane-world model in this work) and obtain the spherically symmetric 
equations in this scenario. For a recent and comprehensive review of the phenomenology of DGP 
cosmology, the reader is referred to [13]. 

In the past two decades, various attempts have been made to understand the nature of the so 
called dark energy and dark matter which would be required to explain many observational data. The 
question of dark matter is presently one of the most intriguing aspects in astrophysics and cosmology. 
There arc some compelling observational evidence for the existence of dark matter for which, the 
galaxy rotation curves and mass discrepancy in cluster of galaxies are two prominent examples [14]. 
Galaxy clusters are the largest virialized structures in the universe, and their mass content is supposed 
to be representative of the universe as a whole. The total mass of galaxy clusters can be determined in 
a variety of ways. The application of the virial theorem to positions and velocities of cluster member 
galaxies is the oldest method of cluster mass determination [15]. More recent methods are based on 
the dynamical analysis of hot x-ray emitting gas [16] and on the gravitational lensing of background 
galaxies [17]. The mass determined from such dynamical means is always found to be in excess of that 
which can be attributed to the visible matter. This is known as the missing mass problem. The 
existence of dark matter was not firmly established until the measurement of the rotational velocity 
of stars and gases orbiting at a distance r from the galactic center could be done with reasonable 
accuracy. Observations show that the rotational velocity increases near the center of a galaxy and 
approaches a nearly constant value with increasing distance from the center. The discrepancy between 
the observed rotational velocity curves and the theoretical prediction from Newtonian mechanics is 
known as the galactic rotation curves problem (figure 2). This discrepancy is explained by 
postulating that every galaxy and cluster of galaxies is embedded in a halo made up of some dark 
matter. To deal with the question of dark matter, a great number of efforts has been concentrated 
on various modifications to the Newtonian gravity or general relativity [18, 19, 20]. There are also 
geometrical approaches to address this problem, namely to use modified Einstein field equations, as 
is done in brane-world models [21, 22, 23, 24] or in modified gravity theories [25, 26]. 

In what follows, we first give a brief review of warped DGP brane-world inspired by the C{R) 
term [27]. Next, we introduce the relativistic collision-less Boltzmann equation from which we deduce 
the generalized virial theorem with the aid of the modified Einstein field equations. The geometric 
mass and density of a cluster is then identified in terms of the observable quantities. Then we derive 
the radial velocity dispersion relation in galaxies. Finally, in the last section, we discuss and conclude 
on our results. 
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Figure 2: Observed rotation velocity curve (solid line B) and the prediction from Newtonian theory (dotted line A). 



2 Einstein equations in warped DGP brane-£(i?) model 

We are going to obtain the relativistic virial theorem for galaxy clusters in warped DGP brane->C(i?) 
model. For this end we need, as a first step, to obtain the gravitational field equations. Consider a 
3-brane S embedded in a five-dimensional bulk M. The effective action of a warped DGP-inspired 
jO{R) scenario is given by 



S = Sbulk + Si 



^brane ) 



where 



^bulk 



I ■ 

Jm 



24 



(1) 



(2) 



Shrane= [ d^X-s/^jCbrane{9a0,4') + f d^X-/^^. 



(3) 



Here, X^, {A = 0,1,2,3,4) are five-dimensional coordinates, while x^, = 0, 1, 2, 3) arc the induced 
four-dimensional coordinates on the brane. = SttGq is the five-dimensional gravitational constant, 

(5) 

TZ and Cm are the five-dimensional scalar curvature and the matter Lagrangian in the bulk, respec- 
tively. K = Vij,n^ is the trace of extrinsic curvature on either side of the brane where 11^ is the unit 
vecror perpendicular to the brane hypersurface S [28]. The last term in equation (3) is known as the 
York-Gibbons-Hawking term which provides a through framework for imposing suitable boundary 
conditions on the field equations. Cbrane{9a/3,'^) is the effective four-dimensional Lagrangian, which 
is given by a generic functional of the brane metric ga/s and matter fields. 
The five-dimensional Einstein field equations are given by 



T^AB - -^TZQab = nl 



Tfl + 5{Y)TAB 



where 



and 



^(5) ^ o <J-^m 
^AB- ^^gAB 



GabC!^\ 



= -2 



(4) 



(5) 



(6) 



We study the case where the induced gravity scenario arises from higher-order corrections on the 
scalar curvature over the brane. The interaction between the bulk gravity and local matter induces 
gravity on the brane through its quantum effects. If we take into account quantum effects of the 
matter fields confined to the brane, the gravitational action on the brane is modified as 



(7) 
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where /x is a mass scale which may correspond to the four-dimensional Planck mass, Af, is the tension 
of the brane and Cm presents the Lagrangian of the matter fields on the brane. We note that for 
C{R) = R, action (1) gives the DGP model if = A(5) = and gives the RSII model if = 0. 
We obtain the gravitational field equations on the brane-world as [4] 



(8) 



where the quadratic correction has the form 



1 



T,? + 



1 

12 



24 



and 



c _ --<(5) „A B c D 



is the projection of the bulk Weyl tensor to the surface orthogonal to n^. 
metric, we can also obtain the surface equation as 



(9) 
(10) 

By the variation over the 
(11) 



where Ki^,^ is the extrinsic curvature and t^i, is the total energy-momentum tensor of the brane. 
Equation (11) is nothing but the Israel junction condition [29]. By taking the divergence of the Israel 
equation, we have 



(12) 



where V^^ is the covariant derivative with respect to g^i,. The right hand side can be evaluated using 
Codacci's equation and the field equation of the bulk, giving 



\7 - -'2T^^K'^n^ 



(13) 



Also, we note that the junction conditions are ambiguous and do not define the four-dimensional 
geometry in the unique way. If one changes the surface action (freedom in the choice of boundary 
terms), the junction condition (11) is also changed. This shows that brane-world is lacking the physical 
principle to predict in unique way the surface term and, hence, the emerging brane cosmology [30]. 
One may prescribe the AdS/CFT correspondence to define the boundary action [31] where curvature 
surface terms correspond to usual countertcrms in dual QFT. 

In order to find the basic field equations on the brane with induced gravity described by C(R) term, 
we have to obtain the energy-momentum tensor of the brane t^i, by definition (6) from Lagrangian 
(7), yielding 



-A{R)gf,u + Tf,^ 



S(i?)G^, + e 



'/if! 



where the functions A{R), T,{R) and Qau are defined as 



MR} = ^ 



dR 



and 



e 



ij.i' 



,dC{R) 
dR 



(dL{R)\ 



(14) 

(15) 
(16) 
(17) 
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Inserting equations (14) into equation (8), we find the effective field equations for the four-dimensional 
metric as [27] 



1 + ^4a{r)^r) 



+ - ^nm^RaG"^ - 4 nrnG^" - Tn 43, - f^., (is) 



where 



(T) 



1 1 1 ^ 2 



(19) 



TT 



— -^Gf^aG'l + —GGfj^i, + -g^ivGapG"''^ - —gfj,„G'^, 



(20) 



and 



24- 



■^jlupa ~ ^ (jjlxvTpcr gupTucr gv<jTfxp) + [^/ii^S'po- + Ti^g^^pgiifj 9ixi>9pcr)\ i 



(21) 



(22) 



43<T = \ {9u.uQp<T - 9u.pQva - 9vaQnp) + [^H^gpa + &{9p.pgva - 9nv9pa)\ , (23) 



3'''V"V„(^) 



with T being the trace of the energy-momentum tensor and is 

e = 9^"®!,, 

and the effective cosmological constant on the brane is given by 

A4{R) 



1 



A(5) + l4AiRf 



(24) 



(25) 



Wc can also recover the standard four-dimensional C{R) gravity [32] from equation (18) in the limit 
A^^ —7- 0. A alternative possibility in recovering the four-dimcnsional C{R) gravity is to take the limit 
0, while keeping the Newtonian gravitational constant 4 finite [4] (see Appendix A). 
Now we are interested in computing the spherically symmetric solutions on the brane. We consider 
an isolated and spherically symmetric cluster being described by a static and spherically symmetric 
metric 



(26) 



we also assume that the matter content of the bulk is just a cosmological constant A^^^ and the matter 
content of the 3-brane universe is considered to be a cosmological constant plus a localized spherically 
symmetric perfect fluid 



iPb + Pb)ViJ.Vu +Pb9iuu- 



Therefore, the gravitational field equations become 

1 



1 + -4a{r)^{r) 



r2 (' 



■1 + rv + 



\4KR)Ph + MR) - l4Hmo + ^Pb 



U{r), 



(27) 
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1 + -4a{r)j:{r) 



4^{R)]Cl^^^ + 4^{R)\l^^^ + 4k:1^^^ - kI^{R)K:1^^^ + 4^1^^^ + ^ [U{r) + 2P{r)] , (28) 



1 + -4a{r)e{r) 



^ (2A' - 2u - ylJr + 2A"r + A'V) = \4^{R)Vh - 2A^{R) + \4A{R)e 



+ 



^ (pg + 2pfePb) - 24Y.{R)K^P + 2/«^S(i?)27r^(^) + 2/«^/C^(®) - 2/«^S(iZ)/C: 



2(G) 
2 



4 2(8) 



+ 2K|7r2^^-^ + 



[t/(r)-P(r)], 



(29) 



where U(r) and P(r) are the dark radiation and dark pressure. Also /C^P, Ifdj^v and /C^®^ are defined 



as 



Now, using these relations and adding the gravitational field equations (27)-(29) we find 



{l + \4A{R)nR)) ^'"{j- 



A' AV a 



A'- 



^k|^[^0(G)_^1(G)_2^ 



2(G) 1 



2^2(0)j_Ac|S(E)2^^O(G)_^l(G) 



^5 ,0(0) 



1(0) 



2.,^(®)] + 



6 



U{r), 



(30) 



where 



^0 -^1 2A:2^ ^ = TT^e [— (1 - e ) + — (1 - e ) \-pbXi^-2pbX - pbX 



12 



4/96 A' 2p6Z/' 6p6Z/' 



(31) 



^O(G) _ ^l(G) _ 2;C2(G) ^ 

U 1 2 r)A 



24 
24' 



^2 /j^3 



-2i/ 



2z/ 



4 4 d-^c 



dR^ 

/\/2 



—R 



TT-R 



o xu. ,,,2 14i^'A' 8A" ^, ,2 22 A' 4A'2 



- 2z.'A" 



■ + 



,'2 ^2 r A uy A,. I 

+ (2A'2 _ 4 + + ^ _ 2AV' + 4A")i?' 

+ ^ _ 2AV + 2A'2 - 4 + 4A")i?''^}, 
r ait'* 



(32) 



^0(0) _ ^1(0) _ 2^2(0) ^ 



d^C d^C 
+ 2(2p6 + 3pb)P"^ + 2(2p6 + 3p6)i?'2_] 



(33) 



6 



and 



TT, 



0(G) 



1(G) 
'1 



1 

12 

2u'\ 



+ 



/,,/2 



12 



r 

4y 4z/ y2 



2 

-4 



4i^' A'2 4A' 2A" 



1^4 1*^3 jj^2 ij-*3 



AV 



+ 



(34) 



0(0) _ 1(0) _ 2 2(0) _ 1 -2. 4 o/y^ 



(2i?" - ly'R' + 



+ 2R 



i/2 



(35) 



eg - e} - 2ei = /x^e- 



2 2 r UR^ dR^ 



(36) 



where a prime represents differentiation with respect to r. In the next section, we will investigate the 
influence of the bulk effects on the dynamics of the galaxies in modified DGP brane-world model. 

3 The virial theorem 



To derive the virial theorem in the context of the model discussed above we have to first introduce 
the following frame of orthonormal vectors [33, 34, 35] 



,(2) 



r sm. 95 



(37) 



where c/^^eiT^ei*^ = r]^"^^^^ 
components is written as 



. The four- velocity v^^ of a typical galaxy with v^^v^ = —1, in tetrad 



,(«) 



a = 0,1,2,3. 



(38) 



Now, we write down the general relativistic Boltzmann equation governing the evolution of the distri- 
bution function of galaxies /b, which are treated as identical and collision-less point particles. This 
relativistic Boltzmann equation in tetrad components is given by 



,{0.)pP 



dfE 



+ 



"(a) Q^p 

/b(x^,i;(")) and 7["^ - .W.P 



T(b)(c)^ ^ 



B 



0, 



(39) 



^p;o"^(6)^(c) ^"^^ distribution function and the Ricci rota- 



where fs 

tion coefficients, respectively. Assuming that the distribution function is only a function of r, the 
relativistic Boltzmann equation becomes 



dfB 
dr 



(vldX 
y2 dr 



- cott* vg 



dv. 



dvr r 

dfB 
dv0 



' dfB ^ dfB 



where 



,(0) = 



,(1) 



dv0 
= 0, 

,(3) . 



(40) 



(41) 



Since our metric is spherically symmetric, the coefficient of cot^ must be zero. Multiplying equation 
(40) by mvrdv where dv = ^dv^dvedv^ is the invariant volume element in the velocity space and m 
is the mass of the galaxy, and integrating over the velocity space and assuming that the distribution 
function vanishes rapidly as the velocities tend to itoo, we find 



dr 



p{vl 



+ 2^ 



dX 
dr 



{vD + {vD - 2{v',) =0, 



(42) 
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where p is the mass density and {v^) represents the average value of v^. Multiplying equation (42) 



by Airr and integrating over the cluster leads to 

rR 



J 

Jo 



47rp {vt) + {vi) + {v') 



r dr — - 



R 



—dr 
or 



This equation can be written as 



2K-- 
2 







R p Q\ 

Anr^p {v'^) + {vl) -Q^dr = 0, 



where the total kinetic energy of the galaxies is defined as 

rR 



K 



/ 

^0 



27rp {vi) + {vi) + {v') 



r dr. 



(43) 



(44) 



(45) 



To obtain the virial theorem in our model we must express the energy-momentum tensor components 
in the terms of the distribution function. This is done according to 



= j fBmVf,Vi,dv, 

which leads to 

Pb = p{vt), Pb = p{vr) = P{'"e) 
Now, using these relations and equation (30) we obtain 



pK). 



(46) 



(47) 



4 ^ 
A 



K 

12 



■i^A{R)p[{v^) + {vD + {vD + {vl)] 



MR) + ^P'iiv^f + {vl? + {vl? + {vl?] - |A(i?)[e° -@\- 201] 



+ [x;^^ - - 2}cf^^] - ^ - - 21c 



2(8)1 



+ 



0(e) 



1(e) 



(48) 



In order to obtain the virial theorem we have also to use some approximations. First, consider that 
the galaxies in the cluster have velocities much smaller than the velocity of light. In other words 



{v'r) « {vD « {vl) « {V^) ^ 1. 



Second, consider inside the galactic clusters the gravitational field is weak and we can use the weak 
gravitational field approximation. Then the term proportional to AV and A'^ in equation (48) may 
be ignored. Third, since we are interested in astrophysical applications at the extra-galactic scale, we 
assume that the deviations from standard general relativity (corresponding to the background value 



^ = 1) are small. Therefore we may represent ^ as ^ = l + eg'{R), where e is a small quantity and 



g'{R) describes the modifications of the geometry due to the presence of the tensor 0^^ [25]. Thus, 
assuming that 1 inside the cluster [21], equations (32), (34) and (35) are zero and equation 

(48) reduces to 



1 



1 



1 + -4Xh ( 1 + ^^f^^^g'iR) ) (1 + ^g'iR)) 



A A 

7 + T 



1 



8 



,„ 2 . , d^£ 



- (i + ^/^9'{R)) {R'o + Ik 



dB? 



(1 + eg'{R)) p [^y'iv-i) + 2A'(t;f) + A"r(^f ) + 3i.'(^;2) 



, 1 4 2 



2 



+ 



6 



■U{r), 



(49) 



where = i?o + ©(A'^, AV). On the other hand, for clusters of galaxies the ratio of the matter 

density and of the brane tension is much smaller than 1, p/\b « 1, so that one can neglect the 
quadratic term in the matter density in above equation. These conditions certainly apply to test 
particles in stable circular motion around galaxies and to the galactic clusters. Thus, we can rewrite 
equation (49) as 



Afe 4 2\ 1 d f 2<9A\ _ kI 



6 



6 



1 + ^4t^' = -^p{r) - A4 + ^Pcurvir) + 4Anr) + 3W(r)] + -^Uir), (50) 



2j.2 Qf y Qj- J 



where 



U{r) = - 



pv_ 
6r ' 



(51) 



gj,2 Qj- y Qff J 



(52) 



Pcurvif) 



2p\ I d ( ^dX 



2Ab Afc I dr \ dr 



(53) 



and 



A — 5 



(54) 



which in the limit of J0{R) = R, we have Pcurvir) = 0- Multiplying equation (50) by and integrating 
from to r yields 

which the total baryonic mass and the geometrical masses of the system are given by 



M(r) = 47r / p{r')r"^dr' , 
Jo 



(56) 



and 



rr 

Mcurv{r) = 47r / Pcurv{r')r"^dr' , 
Jo 

4MDGp{r) = STTKjp'^ l\v{r') + m{r')\r"^dr' , 
Jo 



(57) 



(58) 
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487r , . „ '9 



K^Xb Jo 



KiMRsir) = Z2^ l_ U{r')r 'dr . (59) 



Multiplying equation (55) by ^^^''^ and integrating from to R, we finally obtain the generalized 
virial theorem in modified DGP scenario as 

(l + y «|/^') 2K + W + ^A4l + W^urv + Wdgp + Wrs = 0, (60) 

where 



W 

Stt 



/ — ^dM(r), (61) 
Jo r 



and 



/ r'^dM{r), (62) 
Jo 



are the gravitational potential energy and the moment of inertia of the system, respectively. We can 
also define the gravitational potential energy due to the geometrical masses as 



r Mcurv{r)prdr, (63) 
2 Jo 

Wdgp = -^ MDGp{r)prdr, (64) 
2. Jo 

Wrs = -4 r MRs{r)prdr. (65) 
2 Jo 

By choosing C(R) = R, wc can also obtain the virial theorem in warped DGP brane-world as [23] 

(^1 + y '^s/^') 2K + W + ^Ad + Wdgp + Wrs = 0. (66) 

At this point it is worth nothing that for p = 0, we have M^gp = and Mcurv = so that the virial 
theorem in the modified DGP brane-world is reduced to the virial theorem in the RS brane scenario 
and the gravitational energy modified by Wrs which has its origin in the global bulk effect due to 
the 6^1, term as [21] 

2K + W + \a4I + Wrs = 0. (67) 

In the case = A^^) = 0, the virial theorem in modified DGP brane-world is reduced to the virial 
theorem in DGP-inspired J0{R) gravity scenario 

2K + W + Wcurv + Wdgp = 0, (68) 
for C{R) = R, we can also obtain the virial theorem in DGP brane-world as [22] 

2K + W + Wdgp = 0- (69) 

We note that in equations (68) and (69) Wrs = 0, since in a DGP brane-world with a Minkowski bulk 

space £fj,i, = 0. There is difference between our model with references [21] and [22]. Here, the virial 
theorem modified by Wrs, Wdgp and Wcurv which the first is due to the global bulk effect whereas 
the second term has its origins in the induced gravity on the brane due to quantum correction and 
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the third term arises from higher-order corrections on the scalar curvature over the brane. Also, in 
contrast to the warped DGP model when the five-dimensional contribution is frozen in the modified 
DGP model, the geometrical terms from the four-dimensional C{R) gravity can be used to explain 
the virial theorem mass discrepancy in clusters of galaxies, (see Appendix B) 
Now, we define the following geometrical mass 

Mir) = Mnsir) + MDGp{r) + Mcurv{r), (70) 

so that 

W = -^ M{r)prdr. (71) 
2 Jq 

In order to obtain a relation between the virial mass and the geometrical mass M.{r), we introduce 
the radii Ry , Ri and Tl as [35] 

M2(r) 
Jo -^dM{r) 



n = -^^ (74) 

where Ry is the virial radius and Tl is defined as the geometrical radius of the clusters of galaxies. 
By defining the virial mass as 

and using the relations 

the generalized virial theorem (60) is simplified as 

' + -6^'^ )[-w) + [mJ [-n) ■ ^''^ 

Since the dark matter provides the main mass contribution to clusters, one can ignore the mass 
contribution of the baryonic mass in the clusters and estimate the total mass of the cluster by M.. 
As a result, we conclude that 

Mv{r)^M{T)(^^y'\ (78) 
This shows that the virial mass is proportional to the geometrical mass. 
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4 Estimating the geometrical mass A4{r) 

In clusters of galaxies the most of the baryonie mass is in the form of the intra-cluster gas. The 
following equation provides a reasonably good description of the observational date [36] 

-3|3 

P9ir) = Po{l+'^^ ' , (79) 

where rc is the core radius, and po ^^nd /3 are cluster independent constants. 

The observed x-ray emission from the hot ionized intra-cluster gas is usually interpreted by as- 
suming that the gas is in isothermal equilibrium. Therefore we assume that the gas is in equilibrium 
state having the equation of state Pg{r) = ^^'' pg(r), where ks is Boltzmann constant, Tg is the 
gas temperature, n = 0.61 is the mean atomic weight of the particles in the gas cluster and nip is 
the proton mass. Thus, with the use of the Jeans equation [14], the total mass distribution can be 
obtained as [21, 36] 

Now, substitution of the mass density of the cluster gas in equation (79) gives the total mass profile 
inside the cluster [21, 36] 

24:TTkBTal3 , , 

Mtot r = ^ ^— ^ . 81 

On the other hand, using equation (56)- (58) we can obtain another expression for the total mass 

*-P-^ = 47rrVa(r) + ^T^r^ Pcurv{r) + ^f- [V{r) + m{r)] + -^U{ry. (82) 

Since we have estimated the quantities M(ot(r) and pg(r), the expression for the geometric mass 
density can be readily obtained 

2 ^ 2 ^ 2^2 = /^o h + ^ +Pcurv{r) + ^^[V{r) + 2U{r)] + ^U{r). 83 

Finally, substituting above equation into equation (70), in the limit r » rc considered here, we 
obtain the following geometrical mass 



M{r) ~ 



limpKi " 3(1 - /3) 



r. (84) 



Observations show that the intra cluster gas has a small contribution to the total mass [37, 38, 39], 
thus we can neglect the contribution of the gas to the geometrical mass and rewrite equations (84) as 

, , / 2^TTkBTgl5\ 



To estimate the value of Ai{r), we first note that ksTg ~ 5KeV for most clusters. The virial radius 
of the clusters of galaxies is usually assumed to be r2oo, indicating the radius for which the energy 
density of the cluster becomes P200 = 200pcr, where pcr = 4.6975 x 10~^^ hl^gr / cm^ [36]. Using 
equation (84) we find 

Tcr = 9I.33/3I/2 ' K,'Mpc. (86) 

The total geometrical mass corresponding to this value is 

M{r) = 4.83 X 10^^^^/^ (^) " h^^M^, (87) 
which is consistent with the observational values for the virial mass of clusters [36] . 
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5 Radial velocity dispersion in galactic clusters 

The virial mass can be expressed in terms of the characteristic velocity dispersion ai as 

Mv = ^ajRv. (88) 

By assuming that the velocity distribution in the cluster is isotropic, we have (u^) = = 
3{v^) = 3cr^, with ai and ar are related by Serf = a^. In order to derive the radial velocity dispersion 
(7^ for clusters of galaxies in modified DGP model we start from equation (42) as 



2'^ 



0. 



(89) 



On the other hand, by neglecting the cosmological constant the Einstein field equation (50) becomes 



-I , ■^b 4 2 

1 + y /^s/^ 



J_d_ 
2r2 dr 



,dX 
dr 



$pir) + $Pcurv{r) + 4li' [r{r) + 3Z^(r)] + ^C/(r). (90) 



Integrating, we obtain 



A 



6 



,5A 



dr 



Svr 



(91) 



where Ci is a constant of integration. 

By eliminating A' from equations (89) and (91), it follows that the radial velocity dispersion of 
the galactic cluster satisfies the following differential equation 



KlM{r) 
87rr2 



p{r) 



KlM{r] 
87rr2 



with the general solution given by 



/I I ^fc 4 2\ 2 1 

(1 + y/^sM = - - 



k|M (r) 
87rr2 



p{r)dr 



K\M{r) 
87rr2 



p{r)dr 



p{r) 

1 
P 



^1 ( ^ 



%p{r)dr 



C2 

P ' 



(92) 



(93) 



where C2 is a constant of integration. 

2 



For most clusters > 3 and therefore, in the limit r » r^ the gas density profile (79) is given 
Now, substituting this relation for the gas density profile and equation (85) 



by Pg{r) = po (^) 
into equation (93), for /3 7^ 1 we find 

1 



^4, ,2^ 



[1 + fKtp 

and for ^ = 1 we have 



Po^l 



12(1 -^)(1- 3/3) 



r — 



-3/3 



pnip 



1 



___C2 

(1 + 3^) r PO 



3/3' 



(94) 



(1 + ^4f^') 



PqkI Inr _^ ksTg _^ C3 _^ C2 
8 r jUmp r po 



(95) 



with 4C3 = Ci + As wc noted before for A'^'''^ = A^ = 0, the modified DGP model reduces to the 
DGP model and equation (95) for /3 = 1 reduces to [22] 



Ponj Inr _^ ksTg _^ C3 _^ C2 
8 r pnip r po 



(96) 



In figure 3 we have plotted the radial velocity dispersion for the cluster NGC 5813. The numerical 
values of it are in the ranges P = 0.766, Tc = 25Kpc, ksTg = 0.52^6^, r2oo = 0.87Mpc [36] and 
the radial velocity is about 240fcm/s [40]. As one can see the radial velocity dispersion in modified 
DGP brane- world is compatible with the observed profiles and is the same with the warped DGP 
brane-world model. 
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Figure 3: The radial velocity dispersion in warped DGP brane-£(i?) model for the NGC 5813 cluster with 
Ci = 0.503,1.005 X 10^,2.011 x 10* Mq and C2 = 0.02,0.03,0.04 Ml/Kpc^ for solid, dashed and dot-dashed curves, 
respectively. 

6 Conclusions 



In this paper we have obtained a comprehensive version of the virial theorem within the context 
of the modified DGP brane-world model which can be reduced to the virial theorem in RS, DGP, 
warped DGP brane-world models and C{R) gravity theories. To derive the generalized virial theorem, 
we have used a method based on the collision-less Boltzmann equation. For this end we have also 
considered an isolated, spherically symmetric cluster, situated in a space with metric given by equation 
(26) and described the galaxies, which are treated as identical, collision-less point particles, by a 
distribution function which obeys the general relativistic Boltzmann equation. It has been shown 
that the geometrical mass which is due to the non-local bulk effect and the induced gravity from the 
higher-order corrections of the scalar curvature on the brane provides an effective contribution to the 
gravitational energy, equation (55), which may be used to explain the well-known virial theorem mass 
discrepancy in clusters of galaxies. The virial theorem allows to remove the virial mass discrepancy, 
by showing that the virial mass My is proportional to the geometrical mass M, My ~ M{^)^. 
The geometrical mass can be directly related to the observed virial mass, and for typical clusters of 
galaxies, it is of the same order of magnitude as the virial mass. This shows that at the extra galactic 
scale the geometrical mass plays the role of what is conventionally called the dark matter. 

We have also obtained the radial velocity dispersion profile of clusters and used the observed radial 
velocity dispersion for the cluster NGC 5813 as an example to show that our model can account for 
the velocity dispersion of clusters. Finally, it is worth noting that the advantage of the modified DGP 
model is when the five-dimensional contribution is frozen. In this case the geometrical terms from 
the four-dimensional C{R) gravity can be also used to explain the virial theorem mass discrepancy 
in clusters of galaxies [25] . 
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Appendix A: Field equations in four- dimensional C{R) gravity 

In this appendix we derive the four-dimensional limit of equation (18). Substituting A(i?), S(ii!), 
Q{R) and A4(i?) into equation (18), we obtain 



l + ^-^W^(l + ^(i?5i-£) 



dR 



2Afe dR 



1 4, 



12 



Y.^^dR-^^^Ayl^^dR-^^ 



1 + ffi-^^ n 



2A6 ^ dR ^ 



14 



Ir) 



.dC 



dR 



dC\ 

dR)\ ^"^'•'^'^ 



,d£\ 



dR) 



,(G) 



(97) 



Now we can recover the standard four-dimensional C{R) gravity from above equation in the hmit 
At5 0, while keeping the Newtonian gravitational constant Xb finite [11] 



1 



dC 



dR 



dC\ 



where 



dC 



dR 



1 



neglecting the effective cosmological constant we can also rewrite equation (98) as 

1 



G 



dC) 
dR 



2rp . rp{curv) 

"'4-'- ^IV ^ -'- LLV 



where = 4j and the curvature fluid energy-momentum tensor is defined as 



dC 



dR 



dC 



dC 



dR 



(98) 



(99) 



(100) 



(101) 



which is exactly the four-dimensional field equation in C{R) gravity theories [32]. 



Appendix B: Virial theorem in C{R) gravity 



Our aim in this appendix is to obtain the four-dimensional limit of the virial theorem in modified 
DGP model and to compare its results with pure C{R) gravity theories. From equation (49) the virial 
theorem in modified DGP model is given by 



1 + y /^sM 
1 



J_d_ 

2^2 Qj. 



dr 



^4hp{r) - A4 + l4l^'p + Ik 



dR^ 



p p 3i/' 
1- — + — — 
2Xb Afe r 



-eg\R)-AXbf^^l + ^ + {-]^,i-{r^"£ 



1 d 



2Xb Xb / 2r^dr \ dr J 



>dX\ 



+ 4Ar{r) + 3W(r)] + 4^Uir), 



dR^ 



i^\Xb 



(102) 



in the limit — )• 0, with keeping the Newtonian gravitational constant K4 = \4'^b finite, we obtain 



1 4, 2 1 ^ 



or J 



^ni^bpir) - A4 - eg'{R)^4Xbii^ - £g'{R)^4hp'^^-^ 



1 

12' 



4Xbp,^R'^ + -R'o 



dR'^' 



2d\\ 
dr J 

(103) 



neglecting the effective cosmological constant, A4, we can rewrite the above equation as 



1 d 



,dX 



dr 



(1 + ^9'{R)) ( ) = -^P{r) - -^eg'{R) + Z^'c 



dR^' 



(104) 
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where //^ = 4j. As mentioned before we present 3^ as ^ = l+eg'{R) where e is a small quantity and 

g'{R) describes the modifications of the geometry due to the presence of the tensor Q^u- Therefore, 
we have 

i^¥r i''^) " + 47rGp_.(r), (105) 



where 



AirGpcurv ^ -'iTrGp{r)eg'{R) + 



-sg'iR) - \ + ^rC) 



{l + eg'iR)) \ (106) 



is the geometrical energy density. Equation (105) is exactly equation (17) in reference [25]. Also the 
resulting virial theorem from equation (105) is given by 

2K + W + Wcurv = 0, (107) 

where 

kI fRM{r) 



_^Ar^,Mir), (108) 



w 



Wcarv = / Mcurv{r)pn-dr. (109) 
2 Jo 
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